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Abstract 

We study the dynamics of strongly interacting gauge-theory mat¬ 
ter (modelling quark-gluon plasma) in a boost-invariant setting using 
the AdS/CFT correspondence. Using Fefferman-Graham coordinates 
and with the help of holographic renormalization, we show that perfect 
fluid hydrodynamics emerges at large times as the unique nonsingular 
asymptotic solution of the nonlinear Einstein equations in the bulk. 
The gravity dual can be interpreted as a black hole moving off in the 
fifth dimension. Asymptotic solutions different from perfect fluid be¬ 
haviour can be ruled out by the appearance of curvature singularities 
in the dual bulk geometry. Subasymptotic deviations from perfect 
fluid behaviour remain possible within the same framework. 


1 Introduction 

From the first years of the running of heavy ion collisions at RHIC, evi¬ 
dence have been found that various observables are in good agreement with 
models based on hydrodynamics [I] and with quark-gluon plasma (QGP) 
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in a strongly coupled regime |2]. To a large extent it seems that the QGP 
behaves approximately as a perfect fluid as was first considered in J3|. 

It is a challenge in QCD to derive from first principles the properties of the 
dynamics of a strongly interacting plasma formed in heavy-ion collisions and 
in particular to understand why the perfect fluid hydrodynamic equations 
appear to be relevant. 

Even if the experimental situation is still developing and rather complex, 
it is worth simplifying the problem in order to be able to attack it with 
appropriate theoretical tools. Recently the AdS/CFT correspondence [5| 
emerged as a new approach to study strongly coupled gauge theories. This 
has been largely worked out in the supersymmetric case and in particular for 
the conformal case of J\f = 4 super Yang-Mills theory (SYM). Interestingly 
enough, since the QGP is a deconhned and strongly interacting phase of 
QCD we could expect that results for the nonconfining J\f = 4 theory may 
be relevant. We will make this assumption in our work. 

The AdS/CFT correspondence has already been advocated in theoretical 
studies in the context of heavy-ion collisions HUZllElll. Transport coefficients 
at finite temperature have been calculated using the static black hole dual 
geometry and some generalizations jjG], thermalization has been suggested 
to be described by a black hole formation process |7], and proposals have 
been put forward for the gravity dual description of various processes during 
heavy-ion collisions [2] e.g. cooling as black hole motion in the 5th direction. 

In this paper we focus on the spacetime evolution of the gauge theory 
(4D) energy-momentum tensor, and derive its asymptotic behaviour from 
the solutions of the nonlinear Einstein equations of the gravity dual. 

Imposing the absence of curvature singularities in the gravity dual, we will 
show that, in the boost invariant setting (as in [3]), perfect fluid hydrodynam¬ 
ics emerges from the AdS/CFT solution at large times. The corresponding 
asymptotic solution of the Einstein equations is given by formula (13911 of our 
paper. 

The plan of our paper is as follows. In section 2 we review the Bjorken 
hydrodynamics on the gauge theory side. Then, in section 3, we setup a 
general framework of deriving a gravity dual for a given energy-momentum 
tensor on the boundary, based on the holographic renormalization method. 
In section 4 we derive the large proper-time behaviour of the boost-invariant 
gravity duals by solving analytically the corresponding nonlinear Einstein 
equations in the bulk. In section 5 we arrive at the physical solution by 
requiring the absence of curvature singularities. This constraint selects per- 
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feet fluid hydrodynamics in the 4D gauge theory. We close the paper with 
conclusions and outlook. 


2 Bjorken hydrodynamics 

As is well known, a model of the central rapidity region of heavy-ion reactions 
based on hydrodynamics was pioneered in |3j and involved the assumption 
of boost invariance. In this paper we study the dynamics of strongly inter¬ 
acting gauge-theory matter assuming boost invariance. Let us review now 
the picture which will serve as a basis of our theoretical investigation. 

We will be interested in the spacetime evolution of the energy-momentum 
tensor T flu of the gauge-theory matter. It is convenient to introduce proper¬ 
time (r) and rapidity (y) coordinates in the longitudinal position plane: 

x° = rcoshy x 1 = rsinli?/ . (1) 

In these coordinates the Minkowski metric has the form 

ds 2 = —dr 2 + r 2 dy 2 + dx . (2) 

Assuming for simplicity, y —> —y symmetry and translational and rotational 
symmetry in the transverse plane, the energy-momentum tensor has only 
three nonzero components T tt , T yy and T X2X2 = T X3X3 = T xx , which depend 
only on r. Since we are dealing with a conformal gauge theory, T yu is neces¬ 
sarily traceless 

— T tt + —Tyy + 2 T xx = 0 . (3) 

Energy-momentum conservation D U T^ U = 0 gives a further relation between 
the components: 

T —T Tr + T tt + —rfTyy = 0 (4) 

So using relations ©-(ED, all components of the energy momentum tensor 
can be expressed in terms of a single function /(r): 

( /(r) 0 0 ° \ 

T _ 0 -r 3 f/(r)-r 2 /(r) 0 0 

0 0 f(r) + ±T±f{T) 0 

V 0 0 0 f(r) + ±T£f{T) / 

(5) 
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where the matrix T/„ is expressed in (r,y,xi,x 2 ) coordinates. 

Furthermore the function /(r) is constrained by the positive energy con¬ 
dition which states that for any time-like vector the energy in the frame 
whose timelike axis is should be positive i.e. 

T iW t) l t v > 0 . (6) 

Using t M = (\/ s 2 + t 2 v 2 + 2 w' 2 , v, w , w) we are led to the following restrictions 

f(r) > 0 /'(t) < 0 T f\ T ) > - 4 /( r ) • ( 7 ) 

Note that all the above structure © is purely based on kinematics. The 
dynamics of the gauge theory should pick a specific /(r). A perfect fluid or 
a fluid with nonzero viscosity and/or other transport coefficients will lead to 
different choices of /(r). 

The main aim of this paper is to address the problem of determination of 
the function /(r) from the AdS/CFT correspondence. Let us first describe 
two distinct cases of physical interest: 

Perfect fluid: Bjorken hydrodynamics 

Let us assume that the gauge theory matter behaves like a perfect fluid. This 
means that the energy-momentum tensor has the form 1 

= (E + p)u^u u + pq^ (8) 

where v/ 1 is the (local) 4-velocity of the fluid ( u 2 = —1), E is the energy 
density and p is the pressure. Boost invariant kinematics then forces 2 = 
(1, 0, 0, 0) and the comparison with © leads to 

M = 4 (9) 

r s 

which is the result for the ideal relativistic fluid [3] satisfying E = 3 p. More¬ 
over the entropy per unit rapidity remains constant, while the temperature 
cools down as 

T ~ r” 1/3 . (10) 

1 Note that in ref. 0 the Minkowski metric has (+,—,—,—) signature instead of 
(—, +, +, +) that we use. 

2 In the (r, y, X\, X 2 ) coordinates. 
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Free streaming case 


Let us now consider the free streaming case, where the longitudinal pressure 
vanishes. This property is expected to be valid in the first stages of heavy-ion 
collisions when the QCD coupling is small (see ra for further comments). 
Using © this leads to 

f(r) = - (11) 

r 

where e is a dimensionful constant. 

In the following we will more generally introduce a family of /(r) with 
the large r behaviour of the form 

/(r) ~ t~ s . (12) 

Note that using the energy positivity constraint 0 we are led to consider 3 
0 < s < 4. 


3 Holographic renormalization 

Let us now turn to the AdS/CFT correspondence and describe the dual bulk 
geometry corresponding to a given configuration of the gauge theory energy- 
momentum tensor T^. 

According to the AdS/CFT correspondence, vacuum expectation values 
of (a class of) local operators in the gauge theory can be reconstructed from 
the asymptotics of the dual supergravity fields near the boundary HU. In 
the case of the energy-momentum tensor, the dual field is just the metric. 
The reconstruction of the VEV (T /ti ,) from the near-boundary asymptotics 
of the gravity solution has been first studied in [HUS], and in a systematic 
way in ref. | 13 ]. 

Following [0] we consider general asymptotically AdS metrics in the so- 
called Fefferman- Graham coordinates m 

, 9 q Llu dx tJ 'dx u + dz 2 , . 

ds 2 = ^ --- 13 

z- 

Note that this choice leaves, in general, no remaining diffeomorphism (coordi¬ 
nate) freedom since five conditions on the metric have already been imposed. 

3 An interesting limiting case s = 4 satisfies also the constraint but requires a specific 
treatment which is beyond the scope of our paper. 


5 




One then considers solutions of vacuum Einstein equations 4 with negative 
cosmological constant A = —6 (which corresponds to standard AdS 5 [Ti] l and 
their expansion 5 near the boundary at z = 0: 

9nu = 9$ + z 2 g$ + z*gj£ + z& 9 { $ +... . ( 14 ) 


Here gfw is the physical 4D metric for the gauge theory on the boundary. 
In the following we set it to the flat Minkowski metric g$ = r/ /Jiy . Then 
g$ is found to be zero, while gffl is proportional to the VEV of the energy 
momentum tensor 6 : 

{Thu) = const ■ g$ . (15) 

Hence given some metric in the bulk, i.e. a solution of the supergravity 
equations, the VEV of the gauge theory energy-momentum tensor can be 
directly read off. 

In j!4| one also considers the inverse problem, namely how to construct 
general supergravity solutions of the form m- It turns out that one has 
to give as inputs both gj,°J and g$ in order to generate a solution. Einstein 
equations impose on g$ = (T^ u ) just two consistency constraints : 

(T£) = 0 D v (T^) = 0 (16) 

namely tracelessness (since we are in a conformal theory) and energy momen¬ 
tum conservation' in the gauge theory. Then from these data, the Einstein 
equations allow one to recursively reconstruct in principle all higher terms 
9hu in (HU). This general procedure goes under the name of holographic 
renormalization [Til OE]. 

It is crucial to note that Einstein equations by themselves do not impose 
any further local constraints on g$, or equivalently on (T M „). This raises a 
question about the predictive power of the gauge/gravity correspondence in 
this context. It would seem that a-priori any conserved energy-momentum 
tensor gives a viable gravity background. This would be unacceptable from 

4 Note that there is no energy-momentum tensor in the dual gravity construction. 

5 In (1141) additional logarithmic terms might in principle appear (see (I1J), yet we find 
them absent for the cases (Minkowski metric on the boundary) considered in this paper. 

6 For a generic background metric on the boundary gffl, the formula is more compli¬ 
cated im 

' The covariant derivative here is the one for the gauge theory metric gffl. 
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the gauge theory point of view, since the specific form of (T^) should be 
determined by the gauge theory dynamics. 

Therefore in order to proceed further we need to look for a global condi¬ 
tion which would allow us to determine a physically acceptable solution and 
hence a physical energy-momentum tensor profile selected among all formal 
possibilities. One natural criterion is to require the absence of a ‘naked’ 
singularity in the bulk. 

In this respect we are inspired by the AdS/CFT correspondence for gauge 
theories with Nf ^ 0 flavours na EH nn. There, the embedding of a D7 
brane is constructed with a coordinate y 6, behaving asymptotically 8 for p —> 
oo as 

y 6 = m + 4 + • • • (17) 

P 2 

In Eq. m the leading term is the current quark mass m which is fixed 
(alike to the gauge theory metric gjj°J in our case), while the first subleading 
term c corresponds to the quark condensate (ipip') (alike to the gj$ term 
in ©)• A-priori for fixed m one can construct locally an embedding for 
any condensate c. The requirement that the embedding is nonsingular picks 
|18l TQj the unique physical value of the ( , 0'0) condensate. 

We will see that similar reasoning can be applied for the case of energy- 
momentum tensor. One can construct, order by order, a gravity solution for 
any /(r). The requirement of nonsingularity for the dual geometry will allow 
to pick up the physical /(r). 

Before we proceed to describe boost-invariant geometries let us mention 
two examples where exact solutions of the Einstein equations exist for certain 
non boost-invariant energy-momentum tensors. 


Example I : the static black hole 


Let us consider a static isotropic energy momentum tensor with E = 3p = 
const. The corresponding geometry that we obtain from solving the Einstein 
equations with the boundary condition (HI and a metric of the Fefferman 
form 


ds 2 


— A(z)dt 2 + B(z)dx 2 dz 2 

o + o 


(18) 


s For precise definitions of the variables and the geometrical setting see e.g. m- 
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IS 


ds 2 


(i - */4) 

(1 + z 4 /z$)z 2 


dx 2 dz 2 
dt 2 + (1 + z */4)% + % 


(19) 


and the VEV of the energy momentum tensor can be read off from the 
expansion of the metric d 


(44 oc gj£ = 


( 3/^o 

0 
0 

V 0 


4 0 


0 

i/4 o 
0 
0 


0 \ 

0 

i/4 o 

0 l/z 4 ) 


( 20 ) 


The geometry PI does not look very familiar at first sight. However by 
performing a change of coordinates 


z = 


1 + 


( 21 ) 


we can see that it is exactly the standard AdS static black hole 

2 1 — z 4 jzi 2 dx 2 1 dz 2 

ds 2 = - z^dt 2 + — + - 

z z 1 


z 4 /~z 4 z 2 


( 22 ) 


with zo = zq/\/2. In this way, via the Fefferman-Graham coordinates we 
recover the result of mm 

For later reference let us quote the Hawking temperature (equal to the 
gauge theory temperature) 


and the entropy 



7TZ 0 7TZ 0 


Area 

4G® 




(23) 


(24) 


Example II : the planar shockwave 

The second example of an exact solution using our method is the geometry 
dual to a gauge theory shock-wave (on the boundary). Shockwave solutions 
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have been constructed in AdS spaces [HI HQ whose sources were in the bulk, 
or particles on the boundary Et 

If we introduce light-cone coordinates x = t — y and x + = t + y then 
the background dual to the VEV of the energy-momentum tensor 


T _ = yS(x ) 


(25) 


is (in the Fefferman-Graham coordinates) 


^ 2 — dx dx + + yz 4 6(x )dx 1 + dx\ dz 2 


(26) 


One can check that this is an exact solution of the Einstein equations 9 . 

The metric (HU) represents the gravity background dual to a plane shell 
of matter moving at a speed of light which is an interesting model of e.g. 
an ultrarelativistically boosted large nuclei. This dual background may be 
a good starting point to extend the study of saturation effects to strong 
coupling (as it mimicks closely the setup at the origin of the Colour Glass 
Condensate/JIMWLK picture f23j). 

Ultimately one is interested in the collisions of two such shockwaves ap¬ 
proaching along two light-cone directions which is the setting corresponding 
to heavy-ion reactions. We leave this difficult but interesting problem for 
subsequent work. We will concentrate here on an idealized boost-invariant 
description which would correspond to the description of the central rapidity 
region j3J. 


4 Boost-invariant geometries 

Let us now come to the main issue of this paper namely the study of dual 
geometries in the boost-invariant case. We impose boost invariance, together 
with y —» — y symmetry plus translation and rotation invariance in the trans¬ 
verse plane. The most general form of the bulk metric respecting these sym¬ 
metries in the Fefferman-Graham coordinates reads 

, 2 — e“( r ’ z W 2 + T 2 e b ^ T,z ^dy 2 + dx 2 , dz 2 , . 

_ ds = -12-+ WT • ( 27 ) 

9 The same holds for the more general case where S(x~) is replaced by any function 
of x~. 
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The three coefficient functions a(r,z), b(r,z ) and c(r,z) must start off at 
small z as z 4 according to m-m and ©• In this paper we will restrict 
ourselves to the energy density behaving like 

/M = 4 (28) 

T s 

for 0 < s < 4 and we concentrate on the resulting leading behaviour for 
r —> oo. Let us emphasize that there is a lot of physical content also in the 
subleading behaviour and this problem certainly deserves further study. 
First we solve the Einstein equations 


Rfiv 6 g 0 


(29) 


order by order in z as in (d, starting from and following the holographic 
renormalization procedure. We have implemented the iterative procedure 
using Maple |23, to obtain exact coefficients of the power series expansions 

like 


N 


a{r , z) = y^a n (r)F 


4+2n 


(30) 


n =0 


to some order N. This method calls for comments. 

On the one hand this form is difficult to use in order to analyze possible 
singularities in the bulk since these occur at the edge of the radius of con¬ 
vergence and it is difficult to disentangle unambigously whether the effect 
comes from a finite radius of convergence or is a mark of a genuine curvature 
singularity. 

On the other hand, the knowledge of the power series solution helps us 
to find the large r asymptotics of the exact solutions in an analytic form. 
Namely, by analyzing the structure of the power series solutions (HP, we find 
that after introducing the scaling variable 

z 


V = 


T 


s/4 


(31) 


the exact solutions behave like 


a(r, z) = a(v) + O 


b(r,z) = b(v) + O 


c(r, z) = c(v ) + O 


1 

T # 

1 

r # 

1 

r # 


(32) 
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where we denoted by a positive (here unspecihed) power. 

In order to find a(v), b(v) and c(v) in an analytical form we insert the 
metric (EH) into the Einstein equations © and take the limit r —> oo keeping 
v fixed. We obtain the following set of coupled nonlinear equations: 

v(2a'(v)c'(v) + a'(v)b'(v) + 2b'(v)c'(v)) —6 a!{v) —6 b'(v) — 12c' (v) +vc , (v) 2 = 0 

3 vc\v) 2 + vb\v) 2 + 2 vb"(y) + 4 vc"(v) — 6b' (v) — 12c' (v) + 2 vb'(v)c(v) = 0 

2 vsb"{v) + 2 sb\v) + 8a'(v) — vsa'(v)b'(v) — 8b' (v) + vsb'(v) 2 + 

4 vsc"(v) + 4 sc'(v) — 2 vsa'(v)c'(v) + 2vsc'{y) 2 = 0 . (33) 

Taking a suitable linear combination of these equations and integrating, we 
find that the functions a(v), b(v ), and c(n) satisfy a linear relation 

(4 — 3s)a(v) + (s — 4)b(v) + 2 sc(v) = 0 . (34) 

After nontrivial transformations, the remaining equations may be solved giv¬ 
ing the solution 

a(v) = A(v) — 2 m(v) 

b(v ) = A(v) + (2s — 2 )m(v) 

c(v) = A(v) + (2 — s)m(v ) (35) 

where 

A(n) = ^ (log(l +A(s)u 4 )+ log(l - A(s)n 4 )) (36) 

m(u) = (log(l + A(.s) v A ) - log(l - A(s) u 4 )) (37) 

with 

/3s 2 -8s+ 8 

A(S) = V -24-' (38) 

As a cross-check of this solution we have verified that performing a power 
series expansion of m indeed coincides with the scaling t —> oo limit of the 
exact power series solutions. 

Let us first specialize to the two cases singled out in section 2, especially 
since the perfect fluid case will turn out to be the only one physically relevant. 
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Perfect fluid case 


The perfect fluid corresponds to s = 4/3 in (I28[) . Plugging in this value in 
the above equations leads to the following asymptotic geometry 


ds 2 = — 


^ _ eg 


Z 

1473 


1 + 


eo 
3 ' 


- dr 2 + 


- 4/3 


1 + 

x “ 3 ' 


z 

X73 


{r 2 dy 2 + dxj 


dz 2 

+ —y (39) 


where we reinstated the dimensionful parameter e 0 so that /(r) = e 0 /r^. 

Remarkably enough this geometry is of a form similar to the black hole so¬ 
lution m but with the location of the horizon moving in the bulk according 
to 



From the similarity of the geometry © to the black hole solution © , we 
may qualitatively infer the scaling of the temperature i.e. 

T(r) ~ — ~ t~ 3 (41) 

Zo 

and similarly for the entropy per unit rapidity and transverse area 

S(t ) ~ Area ~ r • —r ~ const (42) 

-o 

in agreement with Bjorken hydrodynamics. 

A word of caution is necessary at this stage. It is not clear whether 
it is possible in the general evolving setting to identify in a precise way 
temperature and entropy of such a geometry (see e.g. discussions in the 
context of dynamical horizons in general relativity El). In addition, in the 
AdS/CFT context, a change of coordinates in the bulk involving both r and 
z will modify which point of the horizon lies ‘above’ which point on the 
boundary thus making the identification of a local temperature and entropy 
density problematic. Nevertheless it is quite probable that some approximate 
notions do exist. 

Finally let us note that our geometry © may be a reliable tool to study 
gauge theory observables which are sensitive to the bulk geometry not far 
away from the boundary. 

In the next section we will indeed show that the solution © is selected 
by a criterion of absence of curvature singularities for large proper-times. 
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Free streaming case 

Inserting s = 1 into the equations (EH) we find that the resulting metric is 
no longer similar to a moving black hole even with some different functional 
form of z 0 (r). Namely one gets 


ds 2 = 


(_ (1 + £)^(1 - + (1 + £)i(l - ^)irV-+ 

z 2 

+(1+^(1-A>^i) dz2 


where v = zf \fr. It is qualitatively different from the perfect fluid case, in 
particular it displays singularities or zeroes at v 4 = a/8 in all coefficients. On 
a more quantitative ground we will now perform an analysis of the curvature 
properties of the whole above class of metrics for 0 < s < 4. 


5 Singularities and curvature 

Looking at the general form of (| 351 ) we see that there is a potential singularity 
for v = A(s) -1 / 4 . However as is often the case in general relativity such a 
singularity may be a purely coordinate singularity as indeed happens in the 
vicinity of the static black hole horizon. In order to unambiguously locate a 
physical singularity we calculate a scalar invariant formed out of the Riemann 
curvature tensor. The simplest one, the Ricci scalar R = g 1 ' 1 ' R llt/ is actually 
by definition equal to —20 for any solution of the Einstein equations (Oil as 
can be directly calculated. Let us then calculate the square of the Riemann 
tensor 

in 2 = RT^R^ (44) 

as a probe of curvature singularities. 

It turns out that we cannot directly reach v = A(s) -1//4 for fixed r since 
we have only an asymptotic solution. Therefore we perform the calculations 
of R and fft 2 by taking r —> oo while keeping v fixed. 

Performing first the calculations for R we find that indeed 

«=-2 ° + o(T). (45) 
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Figure 1: The curvature scalar fH 2 calculated as a function of w — v/A (s) 1 / 4 
for the perfect fluid case s = 4/3 (solid line), s = 4/3 — 0.1 (dotted line) and 
s = 4/3 + 0.2 (dashed line). 


This shows that we can trust the leading asymptotic term and thus that the 
asymptotic approximation is self-consistent. 

The calculation of 91 2 gives a nontrivial result: 


iH 2 = 


(1- A(s)V) 


8 „,32 


10 A (s) v 


A(s)V 4 + 42u 2 VA (s) 4 + 


+ 112 u 24 A (s) 4 - 112 u 24 A (s) 4 s + 36 v 20 s 3 A (s) 2 - 72 v 20 s 2 A (s) 2 + 


+828 A (s) 4 v 16 + 288 u 16 A (s) 2 s - 288 u 16 A (s) 2 - 108 u 16 s 2 A (s) 2 + 
— 136 u 16 s 3 + 27 u 16 s 4 - 320 u 16 s + 160 v 16 + 296u 16 s 2 + 36u 12 s 3 + 


-72u 12 s 2 


A (s) 2 v 8 + 42 uV + 112 v 8 - 112 v 8 s + 10 


, 8„2 


+ 0 -w (46) 


r # 


The striking fact is that in the range 0 < s < 4, the above asymptotic 
expression for 91 2 diverges for v = A(s) 1 ^ 4 for all s apart from the perfect 
fluid case s — | (see e.g. fig. 1.). Remarkably enough, this result requires 
a cancellation of the fourth order pole in front of (PI). We checked this 
analytically by performing a Laurent expansion near that pole. The final 
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result for the perfect fluid case is 


IH 2 


perfect fluid 


8(5 w 16 + 20u> 12 + 174w 8 + 20 w 4 + 5) 

(1 + w 4 ) 4 


(47) 


where w = u/A(|)4 = \fSv. fH 2 reaches just a finite maximum of 94 2 = 112 
at v — 1/\/3, which plays the role of the horizon, starting from the boundary 
value of IK 2 = 40 (see fig. 1.). 

This result means that for asymptotic times and for s ^ | one can reach 
arbitrarily large curvatures in the bulk. This violates our criterion of nonsin¬ 
gular bulk geometry. Thus we should conclude that the asymptotic behaviour 
for large r of the gauge theory energy-momentum tensor should be of the 
perfect fluid type. 

Note that our result does not mean that we have an exact perfect fluid. 
The above analysis was done only in the asymptotic r —> oo regime. We 
showed that gauge theory dynamics rules out such behaviours like streaming 
behaviour which have quite distinct asymptotic /(r) 1/t s with s / f. It 
is quite probable, however, that there would be subleading corrections due 
to e.g. viscosity. In order to detect them one would have to perform a more 
detailed analysis. This certainly deserves further investigation. 


6 Conclusions and outlook 

Let us summarize our main results. 

• We propose a general framework for studying the dynamics of matter 
(plasma) in strongly coupled gauge theory using the AdS/CFT corre¬ 
spondence for the J\f = 4 SYM theory. 

• We use tools related to holographic renormalization for constructing 
dual geometries for given gauge theory energy-momentum tensor pro¬ 
files. We illustrate this method with the static black hole and a planar 
shock wave solution. 

• Further imposing boost-invariant dynamics inspired by the Bjorken hy¬ 
drodynamic picture, we derive the corresponding asymptotic solutions 
of the nonlinear Einstein equations. 
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• Among the family of asymptotic solutions, the only one with bounded 
curvature scalars is the gravity dual of a perfect fluid through its energy- 
momentum tensor profile. 

• This selected nonsingular solution, given by the metric (jUD, is similar 
to a black hole moving off in the 5th dimension as a function of the 
physical proper time. 

Let us add some comments on the specific features of our approach and 
results. In this paper we concentrate on looking for solutions of the full non¬ 
linear Einstein equations. It would be interesting to confront this approach 
with the linearization methods of refs. JT|. In particular viscosity terms are 
expected to appear in the study of subasymptotic terms. Note that the pos¬ 
sibility of black hole formation in the dual geometry has been argued in ref. 

[7j. More specifically, the geometry of a brane moving w.r.t. a black hole 
background has been advocated in ref. |J for the dual description of the 
cooling and expansion of a quark-gluon plasma. In our case we could inter¬ 
pret the solution m as a kind of ‘mirror’ situation in terms of a black hole 
moving off from the AdS boundary. Note however that the precise geometri¬ 
cal identification of the full solution would require further work, in particular 
for the structure near the horizon and for subasymptotic proper times. 

The method and solution presented in this paper raise a lot of stimulating 
questions for future investigation. In particular, one would like to address the 
key problem of connecting, on general grounds, local physical temperature 
and entropy of the gauge theory matter to features of the dual evolving 
gravity solution. 

As a natural outlook it would be interesting to study within the same 
framework possible deviations from perfect fluid behaviour through sub¬ 
asymptotic gravity dual solutions. It would be also interesting to somewhat 
relax boost-invariance and study the possible corresponding modifications 
involving rapidity dependence. In all cases we expect that the condition of 
nonsingularity remains essential to select the proper physical solution. Fi¬ 
nally it would be stimulating to address the initial problem of two colliding 
gauge-theory shock waves using the framework of this paper. 
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